Topographical fingerprints of many-body interference 
in STM junctions on thin insulating films 
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Negative differential conductance (NDC) is a non-linear transport phenomenon ubiquitous in 
molecular nanojunctions. Its physical origin can be the most diverse. In rotationally symmetric 
molecules with orbitally degenerate many-body states it can be ascribed to interference effects. We 
establish in this paper a criterion to identify the interference blocking scenario by correlating the 
spectral and the topographical information achievable in an STM single molecule measurement. 
Simulations of current voltage characteristics as well as constant height and constant current STM 
images for a Cu-Phthalocyanine (CuPc) on a thin insulating film are presented as experimentally 
relevant examples. 

PACS numbers: 85.65.+h, 68.37.Ef, 73.63.-b 



I. INTRODUCTION 

Negative differential conductance (NDC) is a funda- 
mental property of two terminal devices since the discov- 
ery of the first tunnel diodei. The realization of NDC 
within an atomic scale device^— can consequently be re- 
garded as a milestone in the process of miniaturization 
which drives the information technology. 

Scanning tunneling microscopy (STM) experiments 
have played an important role in this research field giving 
several examples of NDC observed with a variety of nano- 
junctions. A number of physical scenarios have been pro- 
posed for the explanation of the experimental findings: 
among others the existence of sharp resonances on both 
electrodes^, the voltage dependent increase in the tun- 
neling barrier heigh^^, the orbital matching between 
molecule and tip 11 ' 12 or even just the symmetry match- 
ing between surface states in the substrate and molecular 
states^. Last but not least vibrational mediated NDC 
has also been observed in single molecule devices^ and 
proposed to test position dependent Franck-Condon fac- 
tors in suspended carbon nanotubea^. 

Recently also interference phenomena in single 
molecule junctions have attracted intense theoretical^— 
and experimental 2 ^— investigations. These junctions 
allow to tackle the fundamental question of the quan- 
tum mechanical nature of the electronic transport at the 
nanoscale and exhibit dramatic modulations of the cur- 
rent desirable for applications. The quest of specific fin- 
gerprints of the electronic interference which go beyond 
the bare current or conductance suppression^ remains, 
though, a crucial issue. We establish in this article a 
criterion to identify the interference blocking scenario by 
correlating the spectral and the topographical informa- 
tion achievable in an STM single molecule measurement. 

In a recent publication we have predicted the occur- 
rence of NDC due to interference blocking 2 ^— in an STM 
single benzene junction on thin insulating filrn^. Ben- 
zene, however is not easily accessible in STM experiments 
and it is not obvious to which extent the findings of [32[ 
apply to larger, experimentally relevant molecules, see 



e.g. (Hp. A major result of this article is an analytical 
expression for the current as a function of the applied bias 
voltage encompassing various transport regimes, see Eq. 
(fTTj) and (H2), which provides both the criteria for the 
occurrence of interference blocking NDC and the inter- 
pretation of its topographical fingerprints. Specifically, 




FIG. 1: (Color online) Artistic view of an STM single molecule 
junction. We show in yellow (light gray) the metallic leads 
(tip and substrate), in red (dark gray) the thin insulating film 
and in green (middle gray) the schematic representation of a 
CuPc. 

interference NDC is expected to occur in molecules which 
exhibit an electron affinity (ionization potential) Eq — E\ 
(E-i — Eq) very close to the work function (f> of the sub- 
strate and, due to their rotational symmetry, have an or- 
bitally degenerate anion (cation) many-body eigenstate 
(here Eq, E±i denote the many-body ground state en- 
ergy of the neutral molecule and of the anion or cation) . 
The necessary decoupling from the substrate, originally 
obtained through a thin insulating layer—, can also be 
achieved by combining two different molecules in a mono- 
layer directly adsorbed on a metal surface^. Recently, 
this set-up has been used to demonstrate position de- 
pendent local gating^, thus suggesting an alternative 
possibility to achieve the mentioned interference condi- 
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tions. Topographical fingerprints to identify the interfer- 
ence blocking scenario are predicted both for the constant 
height and constant current scanning modes. Interfer- 
ence is associated in the first case to a flattening of the 
current map in the molecule region with a corresponding 
loss of the characteristic nodal plane pattern (see Fig. [3]) ; 
in the second case to an enhanced sensitivity of the appar- 
ent height of the molecule to the operating current (see 
Fig. [6]). As we will prove later, both phenomena have 
the same origin i.e., in the interference blocking regime, 
the bottleneck process defining the current pattern is a 
substrate and not a tip tunneling event. The analyti- 
cal results apply to a wide class of molecular junctions. 
In particular, we present simulations concerning a CuPc 
junction. 



II. MODEL 

We describe the STM single molecule junction as a 
system-bath model: 
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where -ff mo i is the Hamiltonian for the isolated molecule, 
in which, to fix the ideas, we distinguish a single particle 
component ffo and a two particle component V, both 
expressed in terms of creation and annihilation operators 
{d a(7l S arJ } for the atomic orbitals i\) a where a indicates 
both the site and the atomic species. ff su b and ffti p 
account for the substrate and the tip, respectively, which 
we assume as reservoirs of non interacting electrons with 
different spatial confinement. Finally fftun describes the 
tunneling coupling between the metallic leads and the 
molecule: 



t^a c , dp 

k£m 



(2) 



where \ = S,T indicates the substrate or the tip, k the 
momentum and a the spin of the electron in the lead. 
Due to their rotational symmetry, the molecular orbitals 
are classified using the projection £ of the angular mo- 
mentum along the principal rotation axis of the molecule. 
A further quantum number m is introduced to account 
for possible degeneracies in the spectrum of the angular 
momentum. Finally, the tunneling amplitudes t 
the form 



t 



Zt m (Ska\lma), 
Stm(Tk z a\£ma), 



(3) 



where si m is the energy eigenvalue of the single particle 
Hamiltonian ffo associated to the state \£ma). For the 
substrate and tip states we assume the model described in 
[32| : a three dimensional momentum is necessary for the 
extended substrate states while only the momentum in 
the transport direction characterizes the tip states which 
are confined in the x and y direction around the tip po- 
sition. 
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FIG. 2: (Color online) Left panel: Current through a CuPc 
single molecule junction as a function of the substrate (and 
tip) work function <f>o and of the sample bias Vb- The tip 
apex position is assumed at (x, y,z — d) — (+5, —5, 7) A with 
the origin taken on the metal-insulator interface and in cor- 
respondance of the center of the molecule, and d being the 
thikness of the insulating layer (see Fig. [TJ. The tip and sub- 
strate resonant lines (respectively with positive and negative 
slopes) divide the parameter space into four regions. T (S) 
indicates a region in which the current is proportional to the 
tip (substrate) tunneling rate. Right panel: Current obtained 
from a cut of the left panel plot corresponding to (j>o = i.leV. 
The numbers on the current-voltage plot refer to the current 
maps of Fig. [3] The current scale is the same for the left and 
right panel. 



Our method of choice to treat the dynamics in the 
regime of weak coupling between system and leads is the 
Liouville equation method. We start from the Liouville 
equation for the total density operator pit) of the whole 
system consisting of the molecule, the tip and the sub- 
strate. We focus on the time evolution of the reduced 
density matrix a — Trs+r{p}, formally obtained by tak- 
ing the trace over the unobserved degrees of freedom of 
the tip and the substrate. A detailed discussion and 
derivation of the equation of motion for the reduced den- 
sity operator of the system can be found e.g. i n 21 ' 36 and 
its adaptation to the STM set-up on thin insulating 
films. For a general discussion about the reduced density 
matrix and related equations of motion see also 37 ' 38 

Let us now consider a molecule deposited on a sub- 
strate having a non degenerate neutral ground state 
\NEo£ = £qS z = 0) and an orbitally degenerate anion 
ground state |N+ 1 E\ I = ±£ x S z = ±1/2). Assume that 
the neutral state minimizes the grand canonical Hamilto- 
nian Hq = H — fig A, where /iq = —<f>a is the equilibrium 
chemical potential of the leads, and (f>o the corresponding 
work function. If E\ — Eq ss fig there is a bias window in 
which the transport characteristics are dominated by a 
dynamics which involves the neutral and anionic ground 
states only. In panel a) of Fig. 2] we give a schematic 
representation of the many body states participating in 
the transport and the associated transition rates where, 
for the sake of simplicity, we neglect the spin degree of 
freedom. According to the general theory presented in 
[321 ]. the corresponding generalized master equation for 
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FIG. 3: Constant height current maps calculated for differ- 
ent bias voltages. The color bar on the left (right) hand side 
corresponds to the maps 1 and 3 (2 and 4). The 5 A long 
white line sets the scale of the images. The numbers in the 
maps refer to the biases indicated in the right panel of Fig. [2] 
The current map in the interference blockade regime (map 4) 
appears flat in the molecule region. The characteristic nodal 
planes pattern appears instead much more pronounced at the 
positive and negative bias resonances (map 1 and 3) and even 
in the Coulomb blockade region (map 2). The tip apex is 
placed at 7A above the molecular plane while the substrate 
biases are, respectively Vbi = 0.1153V, Vb2 = —0.5303V, 
V b3 = -0.7201V, and V b4 = -0.9118V. 



where 

In 

ll, m ,A*w(A£) = -=-E^A/J*^m^-AS), 
k 

(6) 

and we have introduced the notation A£ — £ — t$, At' = 
£' — £q for the variation in angular momenta associated 
to the tunneling process. 

Due to the rotational symmetry of the molecule and 
the different spatial confinement of the leads, the rate 
matrices acquire the form: 

RAi,Ai ,= R $6.1, At' i 

rp rp ( ,Al- Ad \ (7) 

Rai, At> = R ex P ( ~ l ^ 4>Ai I , 

where we did not write for simplicity the energy depen- 
dence of R and the energy and tip position dependence 
of R T and of the phase (j>Ai- Moreover, the latter is de- 
fined as 

PAe = arg \J2 t lAi,n( NE oWAi m r\N+lE 1 £T) . 

. (8) 

Due to their particular structure, the rate matrices ([7]) 
are both diagonalized by the same basis transformation. 



the reduced density matrix in the angular momentum 
basis reads: 
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where t, £' and £" — ±£i , span the angular momenta 
of the anionic ground state and AE = E\ — Eq is the 
energy difference between the anionic and neutral ground 
states. Moreover f£(x) is the Fermi function for the 

lead x, fx( x ) : = f( x ~ Mx) and ; = 1 ~ fx( x )- 

Note that we assume an asymmetric potential drop where 
Ht = A*o — ceV& with c = 0.87 and /is — A*t = eV&. The 
rate R$Le t>-t Q is defined as: 
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FIG. 4: Schematic representation of the many-body states 
participating to the transport. On the vertical axis we re- 
port the grand canonical energies E' := Eq — N/io and 
E[ := E\ — (N + l)/io, being /io the equilibrium chemical 
potential for the leads. In panel a) we adopt the angular mo- 
mentum representation while in panel b) the decoupling basis 
is introduced for the anionic states (see main text for details). 

While the substrate rate matrix is invariant under 
whatever unitary transformation, the tip rate matrix ac- 
quires a peculiar diagonal form since one of its eigenval- 
ues vanishes. The basis transformation, within each spin 
sector of the anionic ground state, reads 



K-7 
Id) 



V2\e 



(9) 



and it depends on the position of the tip via the phase 
4>Ai- Due to the diagonal form of the rate matri- 
ces, in this basis the dynamics is described only by 
means of populations. In particular the decoupled states 
|N + l-Eidr) are only coupled to the neutral ground 
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state \N Eq £q 0) via substrate-molecule tunneling events. 
Both tunneling couplings are still open instead for the 



J 



coupled states |N + 1 E\ c r) , see panel b) in Fig. [4] The 
corresponding master equation reads: 
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(10) 



where for simplicity we have omitted the arguments (A.E) 
of the Fermi functions and the tunneling rates R x and 
suppressed the indexes Eq, £q and E\ in the elements 
of the density matrix. The stationary current flowing 
through the STM junction is calculated as the average 
(I s ) = Tr{cr stat /s} = ~{It) where cr stat is the station- 
ary solution of Eq. (flUt and I x are the current operators 
which are directly obtained from Eq. (Q| following, for 
example; 21 ! 32 . Despite its simplicity, Eq. (| 10[) describes 
the system in a variety of different regimes which leave 
their fingerprints in the current voltage characteristics 
and current maps. 



III. RESULTS 

Given Eq. (flU)) . the stationary current flowing through 
the system is found in closed analytical form, what rep- 
resents one major result of this work. It reads: 



I(Rt ip ,V h ) = 2eR s f+a N 1 



r c N+lT 
N+lr 



(11) 



where e is the (negative) electron charge and 
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Depending the rate R T on the tip position and the bias 
and the Fermi functions on the bias, both topographical 
and spectral information is embedded in Eq. (| 1 1|) . 

In the right panel of Fig. [5] we report the IV charac- 
teristics calculated for a Cu-Phthalocyanine on a metal- 
insulator substrate (a 7A thick insulator with relative di- 
electric constant e r = 5.9) with an effective work function 
(f>s = A.leV . We set up the single particle Hamiltonian 
for the molecule in the tight binding approximation and 
calculate the hopping terms following the Slater-Koster 
scheme 3 ^. Moreover, we adopt the constant interaction 
approximation and assume a charging energy that fits 
the experimentally evaluated electron affinity Eq — E\ of 
CuPc of 4eV. 



At low bias the current is suppressed by Coulomb 
blockade. As the bias increases on the positive side (con- 
ventionally under this condition electrons flow from the 
tip to the substrate) the current undergoes a sudden 
jump corresponding to the opening of the neutral-anion 
transition at the tip- molecule interface (E\ — Eq = /it)- 
On the negative bias side the Coulomb blockade is also 
lifted, but this time at the substrate resonance point 
{E\ — Eq = fJ-s) and the current shows a sharp peak whose 
width scales with the temperature (k^T = QmeV in all 
presented plots). At higher negative biases the current 
is blocked due to interference and the decoupled anionic 
state is the sink of the system. A crucial condition for the 
interference blocking to occur is that E\—Eq <C Eq—E_±, 
ensuring that the substrate-molecule anion resonance an- 
ticipates the tip-molecule cation one which would other- 
wise dominate the transport characteristics. 

Analogous interference blocking involving degenerate 
manybody states has been encountered in a variety of 
systemsiSiSi^i 3 ^. Nevertheless the STM set up de- 
scribed here uniquely allows to correlate the interference 
current blocking with specific topographical fingerprints. 
In Fig. [3] we present different constant height current 
maps (the tip is positioned always 7A above the molecu- 
lar plane) corresponding to the different points labeled in 
the right panel of Fig. [5] Maps 1 and 3 are calculated for 
the tip and substrate resonant tunneling conditions while 
maps 2 and 4 for the Coulomb and interference blockade 
regimes, respectively. Striking is the flattening of the 
current map obtained in the interference case (map 4) if 
compared to all other regimes. Signatures of interference 
can be clearly seen also in the current vs. tip-molecule 
distance represented in Fig. [5] The four traces corre- 
spond to the four different biases conditions indicated 
with the numbers 1 to 4 in the right panel of Fig. [2] and 
the tip is in the same xy position. At large tip-molecule 
distances all traces show the exponentially decaying be- 
haviour typical of the STM measurements (roughly 1 or- 
der of magnitude decay per A). At shorter distances, 
all curves saturates due to the form of the p z orbitals. 
Contrary to the others, though, the curve corresponding 
to the interference blockade regime (case 4) saturates at 
larger distances and shows a wide plateau. For this rea- 
son it even crosses the Coulomb blockade trace (case 2) at 
Az = 7 A, consistently with the result of Fig. [5J Finally, 
we also present in Fig. [S] several constant current topo- 
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FIG. 5: Current vs tip-molecule distance calculated for differ- 
ent biases. The numbers in the legend correspond to the dif- 
ferent cases illustrated in Fig. [3] respectively Vti = 0.1153V, 
V b2 = -0.5303V, Vm = -0.7201V and V b4 = -0.9118V. 
Notice in particular the wide plateau associated to the in- 
terference blockade regime (line 4) and its crossing with the 
Coulomb blockade line for Zu v — d = 7A. 
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FIG. 6: (Color online) Isosurfaces of constant current calcu- 
lated in the proximity of the Coulomb blockade (upper panel, 
Vb — —0.5303 V) and interference blockade (lower panel, 
Vb = —0.9118 V) regimes. The surfaces correspond in both 
cases to the currents: / = 3.15, 3.075, 3.0, 2.925, 2.85 pA. 



graphic maps simulated for different biases and different 
working currents. The surfaces presented in the upper 
panel correspond to the Coulomb blockade regime, while 
the ones in the lower panel to the interference blockade. 
Due to the particular choice of the biases, the apparent 
height of the molecule is exactly the same if we choose 
3 pA as a working current. The shape of the molecule 
is not modified in the interference blockade regime, as 
it is for the constant height current maps (see FigEJ) - 
Yet, striking it is, in this regime, the enhanced sensitiv- 



ity of the apparent height of the molecule with respect 
to the variation of the working current, if compared with 
the same measurement in the Coulomb blockade regime. 
The surfaces presented in Fig. [5] correspond in fact, for 
both cases, to working currents in the range 2.85 - 3.15 
P A. 



IV. DISCUSSION 

All the results presented in the previous section can be 
understood by analyzing the different limits of Eq. (jlll) . 
Let us first consider the Coulomb blockade regime. The 
latter is defined, for Vb < 0, by the inequality E\ — Eq — 
Ms 3> k^T which in turn implies /g <C 1 and <C 1. In 
this limit and under the asymmetry relation R T <C R 
typical of an STM experiment, it is not difficult to prove 
that 



J CB =4ei? T /^ 



ft 
fs 



.4ZL 
fs 



4ei? T /+. (13) 



The current is thus proportional to the tip rate. The 
equality in Eq. (fT3"|) has also a precise physical interpre- 
tation. The charge fluctuations at the substrate lead 
represent the fastest phenomenon {ft / fr ^ 1 due 
to the asymmetric potential drop at tip-molecule and 
substrate-molecule contacts) which sets the ratio between 
the populations of the states to be the thermal average, 
<7 c/d lT / crJV = ft / fs - Finally, the trace sum rule implies: 



a N = I 1 



4ZS. 
fs 



(14) 



The current is determined instead by the slowest pro- 
cess: the tunnelling event |N + 1 E\cr) — > INE^O) to- 
wards the tip. Equation (ITcH) follows due to the presence 
of 2 spin channels and that the tip rate for the coupled 
state is 2R T . Analogously, for Vy, > 0, the Coulomb 
blockade condition reads E\ — Eq — [It S> k^,T and 
the current is again proportional to the tip rate, namely 
/ = —AeR T f^. Thus, the constant height current map 
reproduces the shape of the moelcular orbital encoded in 
R T . 

The interference blockade regime is confined to the neg- 
ative bias and it is defined by the inequality E\ — Eq — 
fis "C —k^T which implies fg « 1 and « 1. Under 
these conditions the current, Eq. (fTTj) . reduces to 



Iib = e 



R fs R fr 
R s fs + R T fr 



(15) 



Equation (fl~5|) tells us, even more clearly if cast into 
the form J^ 1 = {eR s fg)' 1 + {eR T f^)~ 1 , that the cur- 
rent is the result of two competing processes happening 
in series: the thermal unblocking of the decoupled state 
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|N + 1E\At) |N_E o 0) towards the substrate and the 
tip tunnelling event |N + 1E 1 ct) -> |N£; 0). Notice 
that in the system dynamics the two tunnelling events 
are not independent: one cannot happen if the other did 
not happen before. In the interference blocking regime 
fg -C , but, in an STM set up, it typically also holds 
R T <C R s ■ To fix the ideas let us first fix the tip po- 
sition (thus, the ratio R T / R ) and lower the bias, deep 
in the interference blockade, such to fulfill the condition 
R s fg <C R T /y • The current is thus proportional to R s 
and independent of the tip position. This fact explains 
the flattening of the constant height currnet map in Fig. [3] 
and the wide plateau of the current versus tip molecule 
distance in Fig. [5j Nevertheless, as the tip moves far of 
the molecule, the tip rate drops and, as the condition 
R 7 It ^ R S fs ^ s fulfilled, the position dependence of 
the current is recovered (7ib oc R t ). The cross over be- 
tween the two regimes is estimated by the relation: 

R T (R tip , AE) = RSePl&B-vs) (16) 

with (3 — (ksT)^ 1 . For completeness we add that the 
interference blockade is the only regime in which the cur- 
rent looses its canonical dependence on the tip position: 
the current saturates in fact to I = —AeR T for large posi- 
tive biases, it is I — —2eR T at the tip-molecule resonance 
and I = f ei? T at the substrate- molecule resonance. The 
summary of these results and their extension to the work 
function and bias voltage plane is presented in the left 
panel of Fig. [2] where the letters T and S indicate re- 
gions where the current is proportional respectively to 
the tip or substrate rate. 

The enhanced sensitivity of the apparent molecular 
height to the value of the working current of a con- 
stant current scan performed in the interference blocking 
regime can also be explained by analyzing Eqs. (| 13[) and 
(TT5)) . Let us consider a certain value for the working cur- 
rent Iq. Starting from Eqs. (fl3|) and (fl~5|) we can extract 
the equations for the constant current isosurfaces: 

R T = i = K CB (Io,V b ), 

i f i v 1 (17) 

ef T \ eR b f s J 

respectively for the Coulomb blockade and interference 
blockade regimes. If, for a given choice of the parame- 
ters I and V&, it holds K m — Kqb the two associated 
contant current isosurfaces coincide. This is indeed, by 
construction, the case for the bias corresponding to the 
points 2 and 4 in Fig. [2] if the working current is chosen 
exactly as the one in the IV characteristics. Nevertheless, 
for the same choice of the biases, very different sensitiv- 
ity of the constant current isosurface to the value of the 
working current is shown in the interference blockade and 
in the Coulomb blockade cases (compare upper to lower 
panel in Fig. [5]). By analyzing the second equation in 
(TIT)) we can see that Kyb diverges for Iq in the vicinity 



of the interference current eR s fg while Kqb shows a 
completely regular behaviour. As Kyb — > oo the corre- 
sponding isosurface shrinks rapidly as it can be seen in 
Fig. [5] Moreover, the interference current also represents 
in the vicinity of the interference blockade regime, an up- 
per limit for the working current accessible to a constant 
current STM scan. In fact, for Jo > eR s fg the constant 
Kyq turns negative and the second equation in (|17p can 
not be fulfilled for whatever position of the tip. 

Special consideration should be given to the robust- 
ness of the presented effect. Indeed we have presented so 
far the idealized situation in which the rotational sym- 
metry of the CuPc is assumed to be unperturbed with 
a consequently perfect degeneracy of the anion ground 
states. Nevertheless this perfect degeneracy is not a nec- 
essary condition for the occurrence of the many-body in- 
terference effect described in the manuscript. As we have 
already explicitly shown in a previous publication^! the 
interference blocking scenario persists as far as the quasi- 
degeneracy is present i.e. the splitting of the interfering 
energy levels is smaller than the tunneling coupling. In 
fact, if the tunneling coupling is strong enough, the in- 
determination principle does not allow to distinguish be- 
tween the two quasi-degenerate states in the tunneling 
event and interference takes place. Moreover, since the 
tip tunneling coupling is controlled, in an STM experi- 
ment, by the tip position, the interference between quasi- 
degenerate states could be controlled by the tip position. 
The result would be the tuning, with the tip-molecule 
distance, of the negative differential conductance at neg- 
ative bias voltages associated to the interference blocking. 
Finally, for what concerns the effect of the substrate on 
the molecular symmetry we would like to mention that 
strong experimental sensitivity to molecular symmetry 
has been proven for derivatives of CuPc molecules on 
thin insulating films (see Sonnleitner et al. [Io[) suggest- 
ing that an almost complete decoupling of the molecular 
states is indeed a good approximation for these systems. 

Finally, the results presented so far for the CuPc apply 
in general to the class of planar molecules belonging to 
the C nv symmetry group, i.e. invariant under the set of 
rotations of angles z2ir/n, z = 0, . . . , n— 1 around a prin- 
cipal rotation axis perpendicular to the molecular plane 
and to a set of n vertical planes (see Fig.[l|. Their many 
body-states, like the single particle ones, can be classified 
using the projection i of the angular momentum in the 
direction of the main rotational axis (conventionally the 
z axis) that we introduced in Eq. @. The generic many 
body eigenstates of H m can thus be written in the form 
\N E £ S z ), where N is the particle number, E the energy, 
S z and t respectively the projections of the total spin and 
of the angular momentum in the z direction in units of 
fi. The state |N Ei S z ) transforms under a rotation of an 
angle <p = z2n/n around the main rotation axis as: 

R^NEtSz) = e i<p( - e+s ' ) \NE£S z ). (18) 

where R^ is the rotation operator. Consequently it 
is not difficult to prove that £ is an integer number 
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and — j < £ < ^ for C nv molecules with even n and 
— < £ < when n is odd. Since C nv admits at 
maximum bidimensional irreducible representations, we 
conclude that the states with opposite £, connected by 
the reflection operation through the n vertical planes, 
have symmetry protected degeneracy and only states 
with £ = (for even or odd n) or £ — 0, S (for even 
n) are non degenerate. 

V. CONCLUSIONS 

By studying the transport characteristics of an STM 
single molecule junction on a thin insulating film we iden- 
tify in this article a class of molecules that should present 
strong NDC and interference blocking features. More- 
over we establish a criterion to identify the interference 
blocking scenario based on topographical fingerprints. In 
particular, for biases in the vicinity of the interference 
blocking regime, a flattening of the molecular image in 
constant height and an enhanced sensitivity of the appar- 
ent height to the working current in the constant current 
mode are expected. The robustness of the effect is en- 
sured by the observation that quasi-degeneracy and not 
exact degeneracy of the interfering many-body states is 
the necessary condition for the persistence of the phe- 
nomenon. 
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R a = (x a ,y a ,d) T . The corresponding wavefunctions are 
approximated by contracted Gaussian orbitals g2p(f) and 
93d(r) to simplify the calculation of the overlap integrals. 
The definition of the Gaussian orbitals, their contraction 
coefficients d%, e.% and their exponents aj, hi can be found 
i n 41 ' . The orbitals used in this paper then are given by: 

Pz(r) = n 2p f ■ e z g 2p (r), (A3) 

for a p z orbital. A d xz orbital then accordingly reads: 

d xz (r) = n 3d r- e x f- e z g 3 d(r), (A4) 

The parameters n 2p and 713^ are ensuring normalization. 
The electronic states of the tip and the substrate are 
given by |(x = T)ka) or \ (x = S)ka), respectively. Their 
wavefunctions can be expressed in the following form: 

f x (x,y,z) =ip*(x,y)iljX(z), (A5) 

where ip x (x, y) is given by plane waves for \ — S, or by 
the wavefunction of the groundstate of a twodimcnsional 
harmonic oscillator for % = T. The wavefunctions i>\{z) 
are the exponentially decaying parts of the solutions of 
one-dimensional finite potential wells: 

iPl(z) = n s ± e~ KSZ and ifi{z) = n\ e KT(z -^ } , (A6) 

where accounts for normalization and k x is given by: 

K x = \J^f(-4 ~ (A7) 

For the sake of reproduction, the different contributions 
to are listed in the following. 



Appendix A: Tunneling rates and overlap integrals 

The derivation of the tunneling rates, up to a number 
of small differences, is following the example given ir* 3 - 2 ^ 
In their most general form, they are given by: 



r x = 2tt 

im.l'm' t 



E(^J*^w^-A£) (Ala) 



where 



m'cr), (Alb) 

a/3 



M^(AE) = |^^ - AE){ a a\ X ka){ X k<i\Pa). 

k 

(A2) 

The coefficients (£ma\aa) and the energies S£ m are ob- 
tained by diagonalizing the single particle Hamiltonian of 
the molecule, which is set up by using the Slater-Koster 
tight-binding approximation 3 ^. The state \aa) denotes 
an atomic orbital located at site a with position vector 



1. Substrate- molecule tunneling rates 

For two p z orbitals located at sites a and (3, M^p reads: 



s _ 47f4 2 



^-v didj 



S 1 j.S 

■ F +<t>o de 



Jb(fcsfe|)e-^ (o *" 1+ °i 1) 



'0 \[tz 

x F(cii, K S ,~d)F(aj, ks, -d) 



(A8) 



where J n (x) is the n-th order Bessel function, R a p — 
R a - Rp and k s = <J^£(AE -£§ — e z ). The func- 
tion F(a, k, x) results from the overlap of "0* ( z ) with an 
atomic Gaussian orbital and is given by: 



F(a, k, x) 



2a 



— erfc 



K + 2ax 



(A9) 



8 



Here, erfc(a;) is the complementary error function. Con- 
sequently, it follows for a d xz orbital located at R a and 
a p z orbital at Rp : 



2^ 4 
/i 3 



~fc3~ V — n 3dn 2p cos 9 a p ^ p— 



one single integration in energy to evaluate in order to 
obtain the tip-molecule tunneling rates. Because of this, 
they are much more straightforward to calculate than 
their substrate-molecule counterparts: 



Po /-.<-, _ fsl _1 _l s 

d£,-P=Ji(Miia/j|)e-* (6 « } 



x F(bi,K S ,-d)F(a j7 K S ,~d), 



(A10) 



where a/ g is the polar angle of the planar component of 

the vector R a p. From M^ xz j3 one can obtain M% g .p 
by exchanging the corresponding parameters and by re- 
placing cosOa/3 with s'm9 a /3. The expression M^ xz . /3yz 
vanishes exactly due to symmetry reasons, and finally 
M i XZ ;f3,x Z is g iven b Y : 



M T - 2lT 
2tt 



de z D(e z )(aa\Tka)(Tka\/3a)S(eJ - AE) 



J tip 



ft 2 V 2 AE - el 



{aa\Tka)(Tka\/3a). 



(A12) 



a,xz\j3,xz 



2h? 



'■.3(7 



de z ks 



- ]\/fS 

- lvJ -a,yz;j3,yz- 



2. Tip- molecule tunneling rates 



(All) 



Due to the fact that the planar energy component of 
the tip wavefunction is fixed at ey = fku, there is only 



The parameter L t i P stems from the one-dimensional den- 
sity of states of the tip and it is cancelled later on by the 
normalization of the tip wavefunction. Another effect of 
the single integration in energy is that all possible combi- 
nations of p z , d xz and d yz orbitals are surviving. In order 
not to go beyond the constraints of this paper we only 
list the overlap integrals needed to construct the matrices 
MJa . After introducing the following parameters and ab- 



1 ab- 
breviations, v 2 — Ay a 



As 



and finally kt = y jjr(?i<-<J — AE), we are able to give the 
overlap integrals between the different orbitals located at 
R a and the tip wavefunction: 



J 



laaVYka) — —nin^f^Kvy^ — - exp ( — 

z — ' a.,- + v A V 



(a xz a\Tka) = n 3d V2nv 3 V* e J exp 

j W + v ) 

(a yz a\Tka) = n 3d V2irv 3 V* e J exp 



{Axi+Ayi) F( aj ,K T ,d-z tip ) 



(Axl+A y 2 a ) F( aj ,KT,d-Ztt p ) 



{Axi + Ay A a ) F(oj,/ST,d-«ti P ). 



(A13) 
(A14) 
(A15) 
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